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Abstract 

Darwin-Foldy nuclear-size corrections in electronic atoms and nuclear radii 
are discussed from the nuclear-physics perspective. Interpretation of precise 
isotope-shift measurements is formalism dependent, and care must be exer- 
cised in interpreting these results and those obtained from relativistic electron 
scattering from nuclei. We strongly advocate that the entire nuclear-charge 
operator be used in calculating nuclear-size corrections in atoms, rather than 
relegating portions of it to the non-radiative recoil corrections. A prelimi- 
nary examination of the intrinsic deuteron radius obtained from isotope-shift 
measurements suggests the presence of small meson-exchange currents (exotic 
binding contributions of relativistic order) in the nuclear charge operator, which 
contribute approximately ^%. 
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Recent measurements in Garching J!]] and Paris || have greatly improved our 
knowledge of the isotope shift between deuterium and normal hydrogen. Due to their 
much increased precision^], these measurements now rival the traditional relativistic 
electron scattering Q for determining the (nuclear) sizes of these isotopes (and their 
differences). This new level of precision has led to a reexamination of many contribu- 
tions to the level shifts || Q and to the calculation of higher-order QED processes. 
Inevitably, a certain amount of controversy has ensued over the best way to proceed 
and over the proper interpretation of various mechanisms |5|, |6|. Our purpose here is 
to discuss these topics briefly from the nuclear-physics perspective, given that these 
measurements have presented nuclear physics with great opportunities. Nothing that 
we say here is entirely new (indeed, much is very old @, [7], §), but we believe that 
the totality casts considerable light on the interpretation and significance of these 
measurements. 



Specifically, (1) we will (briefly) review the physics from the nuclear-physics per- 
spective. (2) We will discuss the conventions (formalism dependence) attendant to 
introducing nuclear size. Although there is no right or wrong way to do this, there 
are consistent or inconsistent ways to proceed and there are ample opportunities for 
double counting. (3) We will make recommendations for avoiding such problems and 
discuss recent electron-scattering results!! [TO, O from this perspective. (4) We will 



i; 



make a first assessment of the d-p isotope-shift data in terms of "normal" and "exotic 
components of the deuteron structure, even though the latter are not yet entirely well 
defined [JI2"] . A new generation ("second generation") of nuclear potentials |13|, [L4|, [L5| 
gives improved insight into deuteron structure, and this will prove useful in reducing 
theoretical uncertainties. 



Relativistic electron scattering has traditionally been the only successful method 
for measuring the sizes of the lightest nuclei ||. Muonic atoms provided significant 
information on heavier nuclei but until very recently electronic-atom measurements 
lacked the necessary precision. Nuclear physics has been investigated primarily using 
nonrelativistic dynamics, but the increasing precision of electron-scattering data in 
the late 1960's and early 1970's led to a reexamination [[7], || |12[ of the ways that 
relativity can affect a nuclear charge distribution. In order to be as specific as possible, 
we will first discuss various options that have arisen in discussing the simpler and 
better-known proton charge distribution, and then extend the discussion to light 
nuclei. We use natural units (h = c = 1) and the conventions and metric (p 2 = m 2 ) 
of Ref. [IB]. We also remove the proton charge, e p , from all currents. 



For historical reasons (analogy with the electron) the electromagnetic structure 
of the proton was introduced in terms of two form factors (i.e., Lorentz scalars): the 
Dirac form factor, Fi(q 2 ), and the Pauli (anomalous magnetic moment) form factor, 
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i<2(g 2 ). The covariant current (normalized to unit charge) is given by ||T6| 



J x = n(P')(7 A Fi(g 2 ) + ^(l 2 ) ^QuHP) , (1) 

where 7 A and a Xu are Dirac matrices, u(P) and u(P') are Dirac spinors, k p is the 
proton anomalous magnetic moment, M is the nucleon mass, -Fi(O) = ^(0) = 1, and 
q = (P 1 — P) is the momentum transferred (by an electron) to the final nucleon (P') 
from the initial one (P). Because q 2 < for scattering kinematics, it is convenient to 
adopt the SLAC convention, Q 2 = — g 2 > 0, thus avoiding inconvenient minus signs. 



It was soon realized that even though F2 primarily describes magnetic proper- 
ties of the nucleon, it also contributes (in a minor way at small Q 2 ) to the charge 
distribution |]T7J, so the Sachs |T3[] charge and magnetic form factors, Ge and Gm, 
respectively, were introduced: 



Ge 
Gm 



F,m-^F 2 [Q 2 ) 
F^Q 2 ) + n p F 2 {Q 2 ) . 



(2a) 
(26) 



In terms of these form factors, the (laboratory-frame) cross section for (massless) 
electron scattering by protons in first-Born approximation is given by the Rosenbluth 
formula 0, |, ^ 



da 



°"Mott 



A (Q 2 ) + B (Q 2 



+ 1 + 



0^ 

AM 2 



tan 2 (fl/2) 



(3) 



where 6 is the electron scattering angle, ciMott is the cross section for a spinless point 
particle, and 



A (Q 2 ) 
B (Q 2 ) 



GUQ' 



1 + -^ 



G 2 



E ■ 



2M 2 



G 2 M (Q Z 



1 + 



4Af 2 



Q^r 2 

2M 2 



(4a) 
(46) 



Equation (3) applies to elastic electron scattering by an arbitrary nucleus, while Eq. 
(4) applies only to spin-^ systems (such as the proton, 3 He, or 3 H ). The form factors, 
Ge and Gm were proposed long ago|T7|, |20| , |[ as alternatives to Gg and Gm, but 
were never popularly adopted. Equation (3) has been written so that A is a form 
factor associated with the charge distribution, while B is analogously associated with 
the magnetization distribution obtained from the transverse (to q) component of the 
(space) current. This division is most transparently performed in Coulomb gauge [0. 
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Often the square bracket in Eq. (3) is rearranged as [A{Q 2 ) + B(Q 2 ) tan 2 (0/2)], but 
then A is no longer associated solely with the proton charge distribution. 



One has the option of describing the proton's structure in terms of (i 7 !,^), 
(Ge,Gm), or (Ge,Gm)- Only the last option correctly gauges the proton charge 
distribution to order iv/c) 2 (or, equivalently, Q 2 /M 2 ). Factors of r = Q 2 /AM 2 and 
i] = 1 + t are of relativistic origin and also affect the proton mean-square charge 
radius, defined in the Breit frame]?], [L7| as (r 2 ) c h = / d 3 x x 2 p(x), where J A = (p, J). 
Further defining (r 2 )i = — 6 F[(0) and (t 2 )e = — QG' E (0), we obtain from Eq. (2a) 

(r 2 )E=(r 2 ) 1 + ^- 2 , (5a) 

while the charge form factor obtained from Eq. (4a) produces 

(r 2 ) ch = (r 2 ) E + (r 2 ) DF , (56) 
where we have defined (r 2 ) ch = — 6G' E (0) and 

The various mean-square radii, (r 2 ) 1; (r 2 )^, and (r 2 ) c h, differ by amounts of order 
(a? 5 ) ~ 0.044 fm 2 , but are formally identical in the nonrelativistic (large-M) limit. 
Note that (r 2 )^ 2 is often called the proton radius, r p pl| . 



The quantity (3/4M 2 ) in Eq. (5c) is the Darwin- Foldy (DF) term [16|, |22|] and 
is obtained by expanding the 1/rj factor in Eq. (4a). This factor is traditionally 
incorporated into the kinematical factors (along with OMott) and the experimental 
data are then used to determine Ge and Gm- That is, by convention, the Darwin- 
Foldy term is not considered part of the proton structure, even though it affects the 
cross section. 



Nevertheless, to order (1/M 2 ) we can easily expand the A = component of 
Eq. (1) to obtain the true charge density. One finds that the covariant form of u 
(normalized to uu = 1) generates a frame-dependent total charge (obtained by setting 
q — > 0). The reason for this is that the wave function normalization factor 

appropriate for this convention is relegated to the phase space (ie., ^^^ ). If on 

the other hand, we incorporate that factor in J A , the phase space is ( z^Ta ) an d the 



total charge is invariant]^, H]. The invariant form of the charge operator [|1(| [22 



is 



i Q 2 i ri .(^Gm — Ge) /„x 
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where the Darwin- Foldy factor (q 2 /8M) is an explicit part of the charge operator, as is 
the spin-orbit interaction (expressed here in terms of the Pauli spin operator, &) . The 
spin-orbit interaction plays a significant role in the isotopic charge- density differences 
of heavier nuclei 0, Equation (6) for the charge distribution is equivalent (to 
0{1/M 2 )) to using the form factor G E . 

This daunting multiplicity of forms extends to the atomic-physics problem, as well. 
The Barker-Glover |24| calculation of {Za) 4 corrections incorporated the Darwin- 
Foldy part of the charge density as a recoil correction of order (1/M 2 ). This is most 
easily seen by examining the expression that serves as the baseline for defining the 
Lamb-shift energy P5|. Writing 

■ ~ 1/2 

f(» -J) = I 1 + " {ZaY , 2 I (><') 



n-j-\ + ^ ] + \f-{Zay 

then for the state of an electron of mass m e specified by quantum numbers (n,l,j), 
we have to order (Za) 4 /M 2 for the two-body Coulomb problem 

E rJi = me+ M + , [/(n , j) _ 1 ]__^^ [/(n , j) _ 1 ] !+ (|^(^_-2 T )| 1 _,, ], 

(76) 

where /i is the usual reduced mass. This equation can be rewritten as 

E rdj = me+ M +M/ ( I „)_ 1 ]__^^ [/ ( n , j) _ 1 ] 2+ (gM(_l T _-2 T ) +£DF , 

(7c) 

where the contribution of the proton Darwin-Foldy (<5/o) term to the atom's energy is 

Edf ~ w 5io • m 

The standard expression IJ for the leading-order nuclear-finite-size correction to the 
atom's energy is 

= — ^5 — v )ch , (7e) 

and using Eq. (5c) for (r 2 ) c h in Eq. (7e) precisely reproduces Eq. (7d). Consequently, 
the DF term in an atom can be alternatively considered as part of a recoil correction 
of 0(1/M 2 ) (Eq. (7b)) or as the energy shift due to a part of the mean-square radius 
of the nuclear charge distribution (Eq. (7e)). 

Thus, this same Darwin-Foldy term is by convention a recoil correction in atomic 
physics (viz., the Barker-Glover formula, Eq. (7b)) and a kinematical factor in electron 
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scattering (viz., the Rosenbluth formula, Eq. (3)). This is perfectly allowable but 
somewhat confusing, since that term is part of the charge density of the proton in both 
cases. It is unfortunately far too late to change these conventions for the hydrogen 
atom. We do not recommend, however, that they be extended to other nuclei. These 
options were extensively discussed many years ago in the nuclear context and are 
clearly formalism dependent (i.e., a theorist's choice). 

Equation (7b) was originally developed for the proton, but has been applied to 
other nuclei. For the deuteron problem Pachucki and Karshenboim|| have argued 
that the DF term for a pointlike deuteron vanishes, and hence £"df should be dropped 
from Eq. (7c). Khriplovich, Milstein, and Sen'kov|| responded that only the fortu- 
itous choice in Ref. of a particular g-factor for the deuteron caused that term to 
vanish, and in general such a term exists. We agree with Ref.H that this DF term 
should not be included in Eq. (7c), but for different reasons. As we argue below 
(and as noted in Ref. |J), the choice of inclusion or not is formalism dependent, 
although in general the term is not vanishing. Any such term is a part of the nuclear 
charge density (see the discussion below Refs. [§, 0), and contributes a part of the 
mean-square radius of that density. Indeed, as we have seen, whether the proton's DF 
term is a recoil correction or a nuclear-finite-size shift is also formalism dependent, 
although its inclusion in the standard expression (7b) is sanctioned by decades of 
consensus. We strongly advocate that nuclear DF terms be included as part of (r 2 ) c h. 

We examine electron scattering from the deuteron, 3 H, 3 He, and 4 He in turn using 
Eq. (3) 0. This is particularly relevant and topical because of the recent re-analysis of 
the experimental electron-deuteron scattering data by Sick and Trautmann || . Their 
derived radius, (r 2 )^ 2 = 2.128(11) fm, is the rms radius of the complete deuteron 
charge density. This is typical of most nuclear calculations, which work with the 
charge density using the invariant convention (although there are some exceptions). 

The deuteron has Z=l and spin 1, which adds another form factor to the "charge- 
like" form factor, G\, and "magnetic-like" form factor, G 2 : the "quadrupole-like" form 
factor, G3. Various definitions and combinations can be used, and we use the notation 
and definitions of Refs. [[27], p8fl . Because the charge- monopole (the spherical part of 
p) and charge quadrupole (the nonspherical part of p) contributions are incoherent 
(unless the deuteron spin is somehow constrained), the Aq function of Eq. (3) becomes 



A (Q 2 ) = Gl + - 



Q 2 G C 



AM 2 

where for small Q 2 the charge form factor Gc is approximately! 



la) 



G c {Q 2 )^G 1 + ^-Q d , (86) 
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while the quadrupole form factor Gq depends on G , i,G 2 , and G^f^]. The static 
deuteron quadrupole moment is Qd = 0.286 fm 2 . Equation (8b) is equivalent to 
corresponding forms in Refs. [^, || |28], [29|, [HJ. Defining (r 2 ) c h = -6G" C (0) and 



= -6(^(0), one finds 



(r 2 ) ch = (r 2 ) 1 -Q 



d ■ 



(8c) 



Note that (r 2 ) ch is the mean-square charge radius, and not (r 2 )i; — Qd provides a 
Darwin-Foldy-type correction to G\, and is only one part of (r 2 ) c h- Because there are 
alternative form factor definitions for the deuteron, there are corresponding alterna- 
tive size definitions. However, (r 2 ) c h is both unique and physically motivated. 



The 31 
as in Ref . 



5 H and 3 He cases (both having spin - |) mirror the treatment of the proton, 
, where their F C (Q 2 ) is the analogue of Ge in Eq. (2) and F c /r] 1 / 2 is the 
complete charge form factor in the invariant representation. Reference |[Llj| , on the 
other hand, uses a charge operator normalized according to the covariant convention 
and their form factor denoted -F c h(Q 2 ) differs from that of Ref. |]10[ by an additional 
factor of 77 1 / 2 [F c h/r] is the charge form factor if one uses the invariant normalization 
convention] . The mean-square charge radius obtained from Ref. [TU| is therefore given 
by (-6F£(0) + 3/4M 2 ), while from Ref. it is (-6F c ' h (0) + 3/2M 2 ). 



For completeness we also consider the spinless nucleus, 4 He . The form factor and 
the invariant form of the charge operator for a spin-0 nucleus are the same to order 
(f/c) 2 , and there are no DF corrections. We find[||, 1 1| @ Bq - 0) 



^W W) ^ w2) ( 1 + (i))- (9) 

and (r 2 ) c h = — 6Fq(0), which is another attractive property of the invariant form. 

Manifest covariance, which emphasizes form factors, is the traditional way to im- 
plement special relativity, but it is not the only one. Lorentz invariance (at least 
to order (v/c) 2 , which is the limit of our interest here) can be implemented by con- 
structing explicit many-body representations of the Poincare group f| [12|, |30f| . In 
this scheme, no part of the charge density is more fundamental than any other. 
Rather, one works with the complete density, including "boost" effects such as the 
Thomas precession pl|, H]. For these reasons (based on common nuclear practice) 



we strongly recommend the convention that the mean-square radius of the complete 
nuclear charge distribution be used when computing energy shifts. This further im- 
plies that no "Darwin-Foldy" pieces of the mean-square charge radius of a nucleus 
should be incorporated into "recoil" corrections. If the latter is nevertheless done, it 
is imperative that this convention be stated explicitly. 
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Figure 1: Deuteron and proton interactions with external electric field (curly line). 
The nucleons are depicted as double lines, while meson exchanges in deuterium that 
lead to binding or electric currents are shown as shaded double lines connecting the 
proton and neutron. Figure (la) shows the proton, (lb) shows the deuteron graph 
that generates the "matter" radius, while (lc) illustrates meson-exchange currents. 
The graph showing the neutron's finite-size contribution (identical to Fig. (lb) with 
the curly line attached to the neutron) is not shown. 



Whatever conventions are adopted for the proton, consistency within the frame- 
work of nuclear physics (which treats nuclei as composed of nucleons) requires that 
the physics of the deuteron (or any heavier nucleus) incorporate Eq. (6). There will 
be other mechanisms allowed by the presence of additional nucleons, as well. Figure 
(la) shows schematically the interaction of a single proton with an external Coulomb 
field. The solid dot on the double line (the proton) indicates the proton's (finite) 
charge density. An identical interaction occurs in Fig. (lb) on that proton inside the 
deuteron, where again the solid dot indicates the full proton charge distribution in- 
cluding the DF term. We have indicated by shaded vertical bars on left and right the 
strong interactions that bind the proton and neutron together to make a deuteron. 
In addition to the proton interaction, the neutron has a finite size that contributes 
via Eq. (5b) [note that (r 2 )^ vanishes for a system with no net charge]. The exter- 
nal field can attach to the neutron in Fig. (lb) in an identical fashion to the proton 
interaction. As well, the spin-orbit-interaction |7], || (last) term in Eq. (6) generates 
a small relativistic correction (r 2 )so i n the bound deuteron (or any complex nucleus). 
Figure (lc) illustrates a generic contribution of the meson-exchange-current (denoted 
by MEC) type||26||, where the flow of mesons that binds the deuteron generates a 



small contribution of relativistic order to the nuclear charge density ||T2 
Putting everything together, we can write for the deuteron 

(r 2 ) ch = (r 2 ) m + <r 2 )? h + <r 2 >? h + (r 2 ) B , (10a) 



S 



or, equivalently, 

(r 2 ) ch = (r 2 ) pt + (r 2 )"h + (r 2 )? h , (106) 
where the part due to the binding mechanism is given by 

{t 2 )b = (r 2 ) SO + (r 2 ) M EC + • • • , (10c) 
and the "point-nucleon" radius of the deuteron is defined to be 

<r 2 )p t = (r 2 > m + (rV (10d) 

The nucleon mean-square charge radii are given by Eq. (5b) [recall that (r 2 )^ = for 
the neutron case]. In addition, (r 2 ) m is the mean-square "matter" radius, obtained 
directly from the square of the deuteron wave function [(r 2 ) m = J cPr\^d{f)\ 2 (r/2) 2 , 
where r/2 is the distance from the deuteron center-of-mass to the proton]. Equation 
(10) is quite general and applies to an arbitrary nucleus if a factor of N (the number 
of neutrons) multiplies (r 2 )™ h and a factor of Z (the number of protons) multiplies 
( r2 )ch, ( r2 )m, and (r 2 ) p ch . The correction due to nuclear binding mechanisms, (t 2 )b, 
has been written as the sum of spin-orbit contributions from the individual neutrons 
and protons via the last term in Eq. (6) and (potential-dependent) meson-exchange 
currents, plus ... .Its presence makes Eqs. (10) a definition. 

In the traditional interpretation of the isotope shift [1]], one calculates (r 2 )^ — 
(r 2 )^ as the measure of the finite-size difference in the isotope shift, where the first 
(deuteron) term incorporates a proton DF term while the second (proton) term does 
not. This difference then includes a term (?" 2 )df from the proton in the deuteron 
that counterbalances a similar term implicit in the Barker-Glover recoil correction 
for the proton contained in Eq. (7b). This has been done consistently [0]. Thus, the 
proton-size effect (including the DF part) completely cancels in the d—p isotope shift. 
This cancellation must occur on physical grounds (see Fig. (1)), irrespective of the 
fact that in the proton case by convention we choose to call the DF term a "recoil" 
correction, rather than a finite-size term. 

At the level of accuracy of Ref. ||, however, this approach is no longer adequate. 
Each nuclear finite-size effect comes with its own reduced-mass correction (see Eq. 
(7e)). The proton finite-size corrections in the deuterium atom and in the hydrogen 
atom differ by 0.9 kHz in the 2S-1S isotope shift from this effect, although it is very 
tiny for the DF part alone. The finite-size correction should be calculated for each 
isotope with the proper reduced mass before they are subtracted. 



Our final topic is a preliminary analysis of the deuteron charge radius in the 
non-relativistic impulse approximation |2B| (i.e., the "matter" radius). The zero-range 
approximation |32| results from neglecting the d-state wave function and replacing the 
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Table 1: Calculation of the deuteron rms matter radius for a variety of potential 
models listed on the left. The full radius for each potential is shown in the first 
column of numbers, followed by the zero-range approximation for that case, and the 
defect mean-square radius (the difference in the squares of those columns). The 
final column combines the defect with the "experimental" value [[T1J of the zero-range 
approximation (1.9847(18) fm) to obtain a prediction for the full matter radius. 



Potential Model 



(ry/^fm) (r 2 )^(fm) A(r 2 ) (fm 2 ) (r 2 )^ 2 (fm) 



Second-Generation Potentials 



Nijmegen (full-rel) 


1.9632 


1.9811 


-.0705 


1.9669 


Nijmegen (nl-nr) 


1.9659 


1.9831 


-.0681 


1.9675 


Nijmegen (nl-rel) 


1.9666 


1.9839 


-.0683 


1.9675 


Nijmegen (loc-nr) 


1.9671 


1.9843 


-.0680 


1.9675 


Nijmegen (loc-rel) 


1.9675 


1.9847 


-.0680 


1.9675 


Reid Soft Core (93) 


1.9686 


1.9866 


-.0709 


1.9668 


Argonne Vis 


1.9692 


1.9865 


-.0685 


1.9674 


First-Generation 


Potentials 






Reid Soft Core (68) 


1.9569 


1.9683 


-.0446 


1.9735 


Bonn (CS) 


1.9687 


1.9871 


-.0726 


1.9664 


Paris 


1.9714 


1.9890 


-.0695 


1.9672 


de Tourreil-Rouben-Sprung 


1.9751 


1.9926 


-.0694 


1.9672 


Argonne Vu 


1.9816 


2.0005 


-.0754 


1.9657 


Nijmegen (78) 


1.9874 


2.0069 


-.0780 


1.9650 


Super Soft Core (C) 


1.9915 


2.0119 


-.0816 


1.9641 



deuteron reduced s-state wave function by its asymptotic form, Ase'^ 1 ", where (3 is 
the deuteron relativistic wave number and As is the s-wave asymptotic normalization 
constant. This excellent approximation overestimates (r 2 ) 1 / 2 by less than 1%. Table 
1 shows a calculation of (r 2 ) 1//2 for a wide variety of first-generation |34], [35], [36[ [37], [38 



39] , |4Cfj (i.e., older) and second-generation potentials (13], [14], [L5[| (i.e., newer ones that 
fit the nucleon-nucleon scattering data from very well to exceptionally well). The full 
(r 2 ) 1 / 2 is followed by the zero-range result for that potential. The residual, A(r 2 ) = 
(r 2 ) — (r 2 )zR, is next. The residual is small and for our second-generation potentials 
spans the range: -0.0695(15) fm 2 . The zero-range result using the best current 



\2 



values of A s (0.8845(8) fm" 1/2 ) and (3 Q is (r 2 ) ZR = A|/(16/? 3 ) = (1.9847(18) fm) 
which combines with the residual just quoted to give our best theoretical value for 
the root-mean-square (rms) matter radius of the deuteron: 



th 



{r 2 )lL 2 = 1-967(2) fm. (11) 
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Table 2: Experimental and theoretical 2S-1S deuterium-hydrogen isotope shifts in 
kHz. The experimental value is given on the left, followed by the theoretical value 
for point nuclei (with no Darwin-Foldy terms included in either non-radiative recoil 
contribution), the sum of nuclear polarization, nuclear Lamb shift and higher-order 
Coulomb finite-size contributions is next, followed on the right by the leading-order 
nuclear finite-size contribution (including all nuclear Darwin-Foldy terms) adjusted 
to produce agreement with the experimental isotope shift. 



experimental 


point nuclei misc. nuclear nuclear size 


670 994 334(2) 


670 999 503.2 19.2 -5188.4 



This result is our baseline, from which deviations signal "exotic" components of 
the deuteron charge density. We can make our own estimate of this deviation by using 
the current experimental value || of the 1S-2S isotope shift: 670 994 334(2) kHz. We 
also use an updated version of the theoretical analysis presented in Ref. []J, which is 
displayed in Table 2. We use the improved 2* ratio of Ref. || (1836.1526665(40)) 
and the 2* ratio of Ref. @ (1.9990075009(8)). We also use the improved deuteron 
polarizability of Ref. |44j| : the proton polarizability of Ref. |45| cancels in the isotopic 
difference. Higher-order (Za) 5 and [Zaf Coulomb finite-size corrections are ob- 
M 



tained from Ref. 
-0.1140(26) fm 2 



The neutron mean-square charge radius is taken from Ref. [57 



48 1 . Using the deuteron 



All other constants are taken from Ref. 
mean-square charge radii defined by Eq. (10), we obtain the experimental value of 
the deuteron point-nucleon radius 



cxp 



<r 2 > 



1/2 
pt 



1.9753(11) fm. 



(12) 



and 



ex P (r 2 >;r - th(r 2 }^ = 0.008(2) fm, (13) 

where the error in Eq. (12) is obtained by compounding a 1.5 kHz — uncertainty, the 2 
kHz experimental uncertainty, an estimated 4 kHz uncertainty in QED calculations^], 
and an (equivalent) 3.5 kHz uncertainty from the neutron charge radius. These results 
are shown in Table 3. On the scale of these uncertainties the DF terms discussed 
earlier are very large for the 2S-1S transition, approximately 45 kHz/A 2 (A is the 
nucleon number), where roughly 5 kHz changes (r 2 )^ 2 by 0.001 fm. 



The atomic results above can be contrasted with the less precise determination of 
(r 2 )p{ 2 using Eqs. (10) and the electron scattering results of Refs. 



21 



/ r 2\V2 
cxp \i I pt 



1.966(13) fm. 



from which we obtain 



exp\' /pt 



th(r 2 )^ /2 = -0.001(13) fm. 



(14) 
(15) 
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Table 3: Experimental and theoretical deuteron radii. The deuteron matter radius 
corresponding to second-generation nuclear potentials renormalized to the experimen- 
tal zero-range approximation and the experimental point-nucleon charge radius of the 
deuteron are shown in the first two columns, followed by the difference of experimen- 
tal and theoretical results. Relativistic corrections to the mean-square charge radius 
from the electromagnetic spin-orbit interaction and from MEC (assuming minimal 
nonlocality) are listed in the next two columns. The final theoretical estimate of the 
charge radius for pointlike nucleons is listed in the sixth column. No uncertainty is 
given in the final estimate because of consistency problems between the MEC and 
the nuclear potentials. 



t h (r 2 )H 2 (M ex P (r 2 ) p ( 2 (fm) 


Diff. (fm) 


(r 2 ) so (fm 2 ) (r 2 ) M Ec(fm 2 ) 


th (r%( 2 (fm) 


1.967(2) 1.9753(11) 


0.008(2) 


-0.0014 0.0159 


1.971 



At this level of precision, the result (15) is null. Equations (10b) and (12) lead to a 
full deuteron charge radius from the isotope shift of 2.136(5) fm, which is consistent 
with the value of 2.128(11) fm from Ref. @. 

Although the result (13) is effectively nonzero, there is one caveat about its sig- 
nificance. The matter radius derived earlier is not entirely well defined. It was shown 
long ago [T2j that to order (f/c) 2 there are 2 unitary equivalences that arise naturally 



in treating relativistic corrections; these are the (pion) chiral-rotation equivalence 
specified by a parameter, /z, and the quasi-potential equivalence (similar to electro- 
magnetic gauge-dependence) specified by a parameter, v. These parameters modify 
the nuclear potential through nonlocal terms, and also modify the nuclear charge 
operator through meson-exchange currents. Because none of the representations cor- 
responds precisely to a nonrelativistic (i.e., momentum- independent) potential, no 
specification of fi and v is possible without performing a consistent relativistic cal- 
culation (at least to order (v/c) 2 ). Since a unitary transformation cannot change 
observables (and hence the zero-range approximation is unchanged), only the defect 
wave function and the defect mean-square radius ((r 2 ) m — (r 2 }zR,) can be changed and 
both will therefore depend on n and z/, as will (r 2 )MEC- Both ((r 2 ) m + (r 2 )MEc) and 
(?" 2 )ch do not. We can stipulate conditions on the potential that will restrict the pa- 
rameters /i and v. One condition is "minimal nonlocality" , which requires the nuclear 
tensor force to be as local as possible and the entire force to be energy independent. 
This is equivalent to \i = and v = 1/2[12|], and bears a rough correspondence to 
Coulomb gauge in atomic physics. Such a representation is probably the closest to 
(but not quite the same as) using the local potentials that are the norm in nuclear 



physics. This representation for the MEC charge operator is well known [12| and 
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produces 

and together with 
one finds the full radius 



(r 2 ) M EC . 

u=l/2 

2 



= 0.0159 W, (16) 



(r')so = -0.0014 W, (17) 

th (r 2 )H 2 = 1.971 fm, (18) 

which makes up approximately half the difference between the experimental value 
and the baseline estimate predicated on nonrelativistic second-generation potentials: 
(( r2 )pt 2 — ( r2 )m 2 ) given m Table 3. Hopefully, the remaining .004 fm comes from the 
difference between a true relativistic treatment of the deuteron and our nonrelativistic 
one that we have supplemented with (somewhat) ad hoc corrections. Our results for 
(r 2 )e are similar to those of Ref. 



In summary, we have reviewed the various ways that nuclear sizes are incorpo- 
rated into electron scattering and atomic calculations. We strongly recommend the 
convention that complete nuclear charge radii be used in calculating atomic energy 
shifts, rather than radii based on arbitrary form factor definitions. A "baseline" value 
of the deuteron rms radius was calculated using nonrelativistic second-generation po- 
tentials to correct the (excellent) zero-range approximation. A value of the deuteron 
rms radius extracted from the d-p isotope shift is .008(2) fm larger than this baseline 
value, some of which is almost certainly due to meson-exchange currents. A complete 
resolution of the problem caused by this difference awaits a relativistic treatment 



of the deuteron dynamics[50] that is of "second-generation" quality, because we are 



dealing with very small size differences. 
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